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©Finding Every Root of a Broad Class of Real, Continuous 
Functions in a Given Interval 

This robust and reliable algorithm is capable of locating the zeros of a continuous, 
nonlinear function. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


One of the most pervasive needs 
within the Deep Space Network (DSN) 
Metric Prediction Generator (MPG) 
view period event generation is that of 
finding solutions to given occurrence 
conditions. While the general form of an 
equation expresses equivalence between 
its left-hand and right-hand expressions, 
the traditional treatment of the subject 
subtracts the two sides, leaving an ex- 
pression of the form/(jc) = 0. Values of 
the independent variable x satisfying this 
condition are roots, or solutions. Gener- 
ally speaking, there may be no solutions, 
a unique solution, multiple solutions, or 
a continuum of solutions to a given 
equation. 

In particular, all view period events are 
modeled as zero crossings of various met- 
rics; for example, the time at which the el- 
evation of a spacecraft reaches its maxi- 
mum value, as viewed from a Deep Space 
Station (DSS), is found by locating that 
point at which the derivative of the eleva- 


tion function becomes zero. Moreover, 
each event type may have several occur- 
rences within a given time interval of in- 
terest. For example, a spacecraft in a low 
Moon orbit will experience several possi- 
ble occultations per day, each of which 
must be located in time. The MPG is 
charged with finding all specified event 
occurrences that take place within a given 
time interval (or “pass”), without any spe- 
cial clues from operators as to when they 
may occur, for the entire spectrum of mis- 
sions undertaken by the DSN. For each 
event type, the event metric function is a 
known form that can be computed for 
any instant within the interval. 

A method has been created for a 
mathematical root finder to be capable 
of finding all roots of an arbitrary con- 
tinuous function, within a given interval, 
to be subject to very lenient, parameter- 
ized assumptions. One assumption is 
that adjacent roots are separated at least 
by a given amount, xGuard. Any point 


whose function value is less than ef in 
magnitude is considered to be a root, 
and the function values at distances 
xGuard away from a root are larger than 
ef, unless there is another root located in 
this vicinity. A root is considered found 
if, during iteration, two root candidates 
differ by less than a pre-specified ex, and 
the optimum cubic polynomial match- 
ing the function at the end and at two in- 
terval points (that is within a relative 
error fraction £/, at its midpoint) is reli- 
able in indicating whether the function 
has extrema within the interval. The ro- 
bustness of this method depends solely 
on choosing these four parameters that 
control the search. The roots of discon- 
tinuous functions were also found, but 
at degraded performance. 

This work was done by Robert C. Taus- 
worthe of SCT, Inc. and Paul A. Wolgast of 
Caltech for NASA’s Jet Propulsion Laboratory. 
For more information, contact iaoffice @ 
jpl.nasa.gov. NPO-46901 


@ Kalman Orbit Optimized Loop Tracking 

This method has application in military GNSS receivers. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Under certain conditions of low signal 
power and/or high noise, there is insuf- 
ficient signal to noise ratio (SNR) to 
close tracking loops with individual sig- 
nals on orbiting Global Navigation Satel- 
lite System (GNSS) receivers. In addi- 
tion, the processing power available from 
flight computers is not great enough to 
implement a conventional ultra-tight 
coupling tracking loop. This work pro- 
vides a method to track GNSS signals at 
very low SNR without the penalty of re- 
quiring very high processor throughput 
to calculate the loop parameters. 

The Kalman Orbit-Optimized Loop 
(KOOL) tracking approach constitutes a 
filter with a dynamic model and using the 


aggregate of information from all tracked 
GNSS signals to close the tracking loop 
for each signal. For applications where 
there is not a good dynamic model, such 
as very low orbits where atmospheric drag 
models may not be adequate to achieve 
the required accuracy, aiding from an 
IMU (inertial measurement unit) or 
other sensor will be added. The KOOL 
approach is based on research JPL has 
done to allow signal recovery from weak 
and scintillating signals observed during 
the use of GPS signals for limb sounding 
of the Earth’s atmosphere. That ap- 
proach uses the onboard PVT (position, 
velocity, time) solution to generate pre- 
dictions for the range, range rate, and ac- 


celeration of the low-SNR signal. The low- 
SNR signal data are captured by a di- 
rected open loop. KOOL builds on the 
previous open loop tracking by including 
feedback and observable generation 
from the weak-signal channels so that the 
MSR receiver will continue to track and 
provide PVT, range, and Doppler data, 
even when all channels have low SNR. 

The KOOL algorithm will also reduce 
the processor throughput requirements. 
This is enabled because the dynamic 
model of the receiver motion is very 
smooth, so that the full physical orbit 
model can be run at a low rate; for exam- 
ple, every 10 seconds. This contrasts with 
the signal tracking loop requirement for 
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a much less complex set of processor ac- 
tivity at 50 Hz, a 500 times higher rate. 
Coarse benchmarks of PVT filter re- 
quirements for processor throughput, 
and the benchmark tracking loop’s re- 
quirements, indicate KOOL tracking 
will require an order of magnitude less 


throughput, considering both its lower 
rate and greater complexity. 

KOOL tracking high-rate models for 
phase and range at shorter times will be 
generated within the digital logic once 
they are primed with model parameters 
from the PVT filter. The onboard oscil- 


lator must be commensurately stable, re- 
quiring (delta F)/F of about 10 -11 over 
times up to 10 seconds. 

This work was done by Lawrence E. Young 
and Thomas K Meehan of Caltech for NASA’s Jet 
Propulsion Laboratory. For more information, 
contact icuffice@jpl.nasa.gov. NPO-46973 


©Development of Jet Noise Power Spectral Laws 

This model can be used in measuring high-temperature steam pipes, leak noise from high- 
pressure pipes, or any device that generates noise by jet exhaust. 

John H. Glenn Research Center, Cleveland, Ohio 


High-qualityjet noise spectral data meas- 
ured at the Aero-Acoustic Propulsion Lab- 
oratory (AAPL) at NASA Glenn is used to 
develop jet noise scaling laws. A FORTRAN 
algorithm was written that provides de- 
tailed spectral prediction of componentjet 
noise at user-specified conditions. The 
model generates quick estimates of the jet 
mixing noise and the broadband shock-as- 
sociated noise (BBSN) in single-stream, 
axis-symmetric jets within a wide range of 
nozzle operating conditions. 

Shock noise is emitted when super- 
sonic jets exit a nozzle at imperfectly 
expanded conditions. A successful scal- 
ing of the BBSN allows for this noise 
component to be predicted in both 
convergent and convergent-divergent 
nozzles. 

Configurations considered in this 
study consisted of convergent and con- 
vergent-divergent nozzles. Velocity ex- 
ponents for the jet mixing noise were 


evaluated as a function of observer 
angle and jet temperature. Similar in- 
tensity laws were developed for the 
broadband shock-associated noise in su- 
personic jets. 

A computer program called “sjet” was 
developed that provides a quick estimate 
of component noise in single-stream jets 
at a wide range of operating conditions. 
A number of features have been incor- 
porated into the data bank and subse- 
quent scaling in order to improve jet 
noise predictions. Measurements have 
been converted to a lossless format. Set 
points have been carefully selected to 
minimize the instability-related noise at 
small aft angles. Regression parameters 
have been scrutinized for error bounds 
at each angle. Screech-related amplifica- 
tion noise has been kept to a minimum 
to ensure that the velocity exponents for 
the jet mixing noise remain free of am- 
plifications. A shock-noise-intensity scal- 


ing has been developed independent of 
the nozzle design point. 

The computer program provides de- 
tailed narrow-band spectral predictions 
for component noise (mixing noise and 
shock associated noise), as well as the 
total noise. Although the methodology 
is confined to single streams, efforts are 
underway to generate a data bank and 
algorithm applicable to dual-stream jets. 
Shock-associated noise in high-powered 
jets such as military aircraft can benefit 
from these predictions. 

This work was done by Abbas Khavaran 
and James Bridges of Glenn Research Center. 
Further information is contained in a TSP 
(seepage 1 ). 

Inquiries concerning rights for the commer- 
cial use of this invention should be addressed 
to NASA Glenn Research Center, Innovative 
Partnerships Office, Attn: Steven Fedor, Mail 
Stop 4-8, 21000 Brookpark Road, Cleve- 
land, Ohio 44135. Refer to LEW-18600-1. 
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